Coupling the inflaton to an expanding aether 

William Donnelljjj and Ted Jacobsorjj 
Center for Fundamental Physics, Department of Physics 
University of Maryland, College Park, MD 20742-4111 

We consider a Lorentz-violating theory of inflation consisting of Einstein-aether theory with a 
scalar inflaton coupled bilinearly to the expansion of the aether. We determine the conditions for 
linearized stability, positive energy and vanishing of preferred-frame post-Newtonian parameters, 
and find that all these conditions can be met. In homogeneous and isotropic cosmology, the inflaton- 
aether expansion coupling leads to a driving force on the inflaton that is proportional to the Hubble 
parameter. This force affects the slow-roll dynamics, but still allows for a natural end to inflation. 
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I. INTRODUCTION 

The space-time continuum may be only a low- 
resolution approximation to a more structured, discrete 
plenum, with a UV cutoff on degrees of freedom. This 
structure might be revealed in high energy collisions or 
subtle effects of long distance propagation, or discrete 
symmetry violations. But it might also manifest itself at 
a macroscopic level as a result of the expansion of the 
Universe. The expansion would presumably require that 
new "atoms of space-time" together with new degrees of 
freedom are created as time goes on. The aim of this 
paper is to explore possible consequences of this scenario 
using an effective field theory description. 

A UV cutoff is inconsistent with locality and Lorentz 
symmetry, so at least one of these properties would have 
to fail if such a cutoff exists. In this paper we study a 
model theory in which locality is preserved, but Lorentz 
symmetry fails. To accommodate this scenario in an ef- 
fective field theory that preserves general covariance and 
the successes of general relativity, one needs to incor- 
porate a Lorentz-violating dynamical field into the the- 
ory. Einstein-aether theory provides a straightforward 
approach to doing this. In that theory, the local struc- 
ture of space-time is described by a metric tensor, as in 
general relativity, but also by a dynamical unit timclikc 
vector field u a , called the aether. The aether defines a 
preferred rest frame at each point of space-time, but pre- 
serves rotational symmetry in that frame. This structure 
would suffice to accommodate a local UV cutoff. For 
a Lorentz violating theory like this to be viable, there 
must ultimately be a reason why conspicuous Lorentz vi- 
olation does not infect the low-energy matter action in 
fiat space-time. For the purposes of this paper, we will 
simply assume that such a reason exists. 

Just as the Einstein-Hubert action is (besides a cosmo- 
logical constant term) the unique lowest order covariant 
term in a derivative expansion of the action, the action 



for "pure" Einstein-aether theory consists, in addition, of 
the four independent two-derivative terms involving the 
metric and the aether. Direct coupling of the aether to 
matter would entail local Lorentz violation and is severely 
constrained by observational bounds [lj. 

Our primary goal is to use Einstein-aether theory to 
model possible consequences associated with the growth 
of the vacuum in cosmology. In particular, we study the 
leading order coupling of a scalar inflaton field tp to the 
local expansion rate. Such a term cannot appear in a 
purely metric theory because the local expansion rate, 
i.e. the Hubble parameter H , cannot be constructed from 
the metric in a covariant way. However, given the aether 
field u a , the local expansion 9 relative to the preferred 
frame is a space-time scalar, 



9 = V a u a . 



(1) 



(In a homogeneous and isotropic cosmology, 9 is related 
to the Hubble parameter by 9 — 3H.) Any scalar oper- 
ator might couple to 9, but here we shall consider just 
a (still mysterious) scalar inflaton, since that would pre- 
sumably dominate during inflation. This of course leaves 
open the possibility that scalar operators composed of 
matter fields could couple to 9. At late times 9 is ev- 
erywhere small relative to particle physics scales, so 9 
couplings would not lead to any observable Lorentz vio- 
lation in particle physics. 

We shall consider a generic coupling of the inflaton to 
the expansion, but the lowest dimension (and therefore 
presumably dominant) coupling takes the form 



dtp = —u a V a tp + total derivative. 



(2) 
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This peculiar leading order term has not, as far as we 
know, been considered previously. It apparently entails 
a violation of time-reversal symmetry, although if the 
rest of the inflaton action is even in tp then time-reversal 
symmetry is preserved when accompanied with tp —¥ — tp. 
We first examine the impact of this new coupling on all 
of the noncosmological theoretical and observational con- 
straints on the theory, and then we study the modified 
inflationary dynamics. A key consequence of this term 
is that the expansion 9 acts as an external force, which 
can either slow down or accelerate the evolution of the 



inflaton. Its effect on the primordial fluctuations has not 
yet been fully investigated. 

Although the particular coupling we examine may be 
new, much work has already been done on cosmology in 
Einstein-aether theory. Inflation with a preferred time- 
like vector was first considered by Gasperini [2]. More 
recently, authors have computed the effect of the aether 
on the spectrum of primordial perturbations [3-6 J , with 
no direct coupling between the aether and the inflaton. 
A certain kind of aether-inflaton coupling has been con- 
sidered in which the kinetic terms in the aether action are 
multiplied by functions of the scalar field [7J [5] • In this 
model the effective gravitational constant depends on the 
inflaton, leading to the possibility of a repulsive phase 
that significantly changes the dynamics of inflation. A 
similar model with a nonnormalized vector field playing 
the role of the aether and scalar fields has been proposed 
as a solution of the cosmological constant problem by 
a similar mechanism 9J. Cosmology with a Lagrangian 
that is an arbitrary function of the usual Einstein-aether 
Lagrangian has also been considered as a possible source 
of dark energy [TUl HJ . 

This paper is organized as follows. In Sec. [TT] the the- 
ory is defined and the equations of motion are given. In 
Sees. Ill IV and [V] we work out the properties of lin- 
earized spin-0 perturbations and the conditions for their 
stability, absence of vacuum Cherenkov radiation, and 
energy positivity. This is necessary in order to estab- 
lish the conditions for viability of the model with aether- 
scalar coupling. In Sec. |VI| these conditions are combined 
with those for the spin-1 and spin-2 perturbations of pure 
Einstein-aether theory, the effect of the new coupling on 
the parametrized post-Newtonian (PPN) expansion is de- 
duced, and some additional observational constraints are 
discussed. In Sec. |VII| the homogeneous, isotropic cosmo- 
logical dynamics is investigated, and we conclude with a 
summary and some further remarks in Sec. |VIII| 



II. ACTION AND EQUATIONS OF MOTION 

Einstein-aether theory is a theory of a metric g ab of 
signature (+, — , — , — ) and a unit timelike vector field u a 
|12j . This theory is defined by the Lagrangian 



-M 2 



(R + K ab cd V a u c V b u d + \{u a u a - 1)) , (3) 



where M = (87rG)~2 is the reduced Planck mass. In 
general relativity, the sign of G can be fixed by requiring 
that gravitational waves carry positive energy density. 
The same turns out to be true in Einstein-aether theory, 
hence we will assume from the outset that G > 0, or 
equivalently M 2 > 0. The three terms in the action are 
the Ricci scalar for g ab , the Einstein- aether kinetic term, 
and a Lagrange multiplier term that forces u a to be a 
timclike unit vector. The tensor K ab cd is given by 



where ci.2,3,4 are the dimensionless free parameters of 
Einstein-aether theory. 

The theory defined by the Lagrangian pi) has a New- 
tonian limit with a modified Newton's constant 



Gjv — k nG, 



k n 



C14 



(5) 



The post-Newtonian expansion will be discussed in 



Sec. VI B The Friedmann equation for spatially flat ho- 
mogeneous cosmologies is also the same as in general rel- 
ativity, but with a renormalized Newton's constant Gq 
in place of G PSE], 



G 



c 



K C G, 



KC 



2 + cia + 3c 2 



(6) 



The spatial curvature term in the Friedmann equation is 
also renormalized. 

We introduce to the theory a scalar inflaton p with a 
potential that depends both on ip and the expansion of 
the aether, 



L v = -V a <pV a <p- 



V(6,<p). 



(7) 



For special cases of the potential V(9, <p) the cosmolog- 
ical dynamics of our model overlaps with other models 
of Lorentz-violating inflation. To see this, we note that 
in a Friedmann-Lemaitre-Robertson- Walker metric with 
the aether aligned to the cosmological rest frame, 



V a M6 = H(g ab - u a u b ), 



(8) 



where H = a/a is the Hubble parameter. It follows 
that the aether kinetic term in pi is proportional to 
9 2 . In the model of Kanno and Soda [TJ, 01,2,3,4 & re 
allowed to depend on a scalar field p. This will have 
the same cosmological dynamics as our model with a 
potential of the form V(9, <p) — f(ip)9 2 . The model of 
Zlosnik et. al. [10 does not have a scalar field, but its 
Lagrangian contains an arbitrary function of the aether 
kinetic term. This cosmology corresponds to a potential 
of the form V(9, <p) — f(9 2 ). This equivalence holds only 
for the background dynamics; the linearized perturba- 
tions of these models will generally differ. 
The equation of motion for ip is 



Up + v v = 0, 



(9) 



where V v denotes the partial derivative dV/d<p. The 
equation of motion for u is 

V b K\ = \u a + c A (u c V c u b )V a u b - Ar 2 V a V 9l (10) 



K ab cd V a u c . The component of Eq. (10) 



where K c 

parallel to u a determines A. Finally we have the Einstein 

equation: 



K ab cd = c l9 ab g cd + c 2 5 a c 5 b d + c 3 S2S b + c A u a u b g cd , (4) 



Gab — S ab + M T ab , 



(11) 



where S ab is the energy-momentum tensor for the aether 
in the absence of the scalar field [15] 



&ab 



'c ( K (a U b) ~ RC (a U b) - K {ab) U C 

ci(V c M a V c u & - V a u c V fc u c ) 

c 4 (u c V c u a )(u d V d u b ) 

(u c V d K d c - c 4 (u c V c u e )(u d V d u e )) u a u b 

lg ab (K c d V c u d ), 



(12) 



round brackets denoting symmetrization. The matter en- 
ergy momentum tensor T ab is obtained by varying the 
matter action with respect to the metric. Accounting 
for the contribution of the Lagrange multiplier term the 
energy-momentum tensor can be written in terms of the 
Lagrangian as 



T, 



5L 



(lb 



&g 



ab 



u -^— c UaUb 



Lg a b 



For the scalar Lagrangian Q this formula gives 

T a b = v a ^v fc ^ - (iv a <£>v°y - v + evg) 9ab 

+ (u c V c V e )(u a u b - g ab ). 



(13) 



(14) 



Note that terms in the action linear in 6 do not con- 
tribute to V — 9Vg, hence they contribute only to the 
part of the stress tensor proportional to the spatial met- 
ric u a u b — g ab . This can be traced to the contributions 
from the Lagrange multiplier term in the action, but 
it can also be understood more directly from the fact 
that the metric enters such terms via the combination 
\J—gQ = d a (y/—gu a ). In order to maintain the unit norm 
condition g ab u a u b = 1, a metric variation must be accom- 
panied by an aether variation 8\\u a = —\{5g m n u m u n )u a . 
The combined 5g ab and 8»u a variations yield 



<5(\/ 3 5" a ) = hV-gu a {-g mn +u m u n )5g r ' 



(15) 



The corresponding contribution to the stress tensor is 
thus proportional to the spatial metric, and corresponds 
to an isotropic pressure in the aether frame. 



III. LINEARIZED PERTURBATIONS AND 
DISPERSION RELATIONS 

We consider linearized perturbations about the flat 
space-time solution with a constant aether field 



gab = Vab, 



<5 °, <P = 0. 



(16) 



To ensure that ( 16 ) is a solution of the equations of mo- 



tion, the potential must satisfy 

v(o,o) = v^(o,o) = o. 



(17) 



In order to study linear perturbations, the potential is 
expanded to second order in the fields 



V(6,<p)n \m 2 ip 2 + ixMOy 



(18) 



where m is the mass of the scalar field. The constant 
\x has dimensions of mass, and we include a factor of 
M to compensate for the canonical normalization of the 
aether perturbations. The expansion ( 18 ) does not in- 



clude a term linear in 9, because such a term is a to- 
tal divergence and therefore does not affect the classical 
equations of motion. The term proportional to 2 is also 
omitted, because it can be absorbed into the c 2 term in 
the Einstein-aether Lagrangian. 



Because the background solution ( 16 ) is invariant un 



der spatial rotations, we can decompose the perturba- 
tions into irreducible representations of SO(3). The prop- 
agating degrees of freedom consist of a spin- 2 graviton, a 
spin-1 aether-metric wave, a spin-0 aether-metric wave, 
and a spin-0 inflaton perturbation. At the linearized 
level, the presence of the scalar field affects only the spin- 
modes. This is because the term 8<p in the action de- 
pends only on the timelike and longitudinal parts of the 
aether field, both of which have spin-0. 

We will adopt the variables and gauge choice of Foster 
|16j . In this parametrization, the seven spin-0 degrees 
of freedom are labeled h 0Q , 7, <f>, /, w , v, tp, and a general 
spin-0 perturbation of the metric and aether can be writ- 
ten as 



gab 



1 + h 00 
d tl -1 

[l + w°,8iv]. 



didji 



da 



didjf) 



(19) 
(20) 



Here labels the time coordinate, latin indices i,j,... — 
1, 2, 3 are spatial, and A = S^didj. We can simplify the 
equations of motion by using the gauge freedom to set 
^ = 7 = [16] . Additionally, we can eliminate w° using 
the constraint u a u a = 1, giving 



w 



nhoo- 



(21) 



With these choices, and taking perturbations to have the 
space-time dependence e'C™ - ^') in terms of the flat space 
coordinates (t, X, y, z), the metric and aether take the 
form 



g a b = diag (1 + h 00 , -1 - k 



(1-i/ioo, 0,0,0). 



1 



1 1.2 

2 



k 2 f,-l 



\k 2 f) , 

(22) 

(23) 



In terms of these fields the linearized versions of the 
scalar field equation, the (0,0) Einstein equation, the 
(1,1) Einstein equation and the sum of the (2,2) and (3,3) 
Einstein equations are 

= -uj 2 tp + k 2 (p + m 2 (p - \k 2 iujnM(f + 4>), (24) 

= -k 2 .f - c 14 h 00 , (25) 

= (1 + c 2 )k 2 uj 2 f + c 123 k 2 uj 2 <f> + 2[iiuj(p/M, (26) 
= (1 + cis + 2c 2 )k 2 uj 2 f 

- k 2 (k 2 f + 2/iqo - 2(1 + c 2 )uj 2 <J)) + i/jLiuip/M. (27) 



where we have adopted the notation C14 = C\ + C4, C123 = 
c\ + c 2 + c 3 , etc. 

In the absence of the scalar field the linearized modes 
all satisfy a linear dispersion relation of the form uj = s.jfc 
where Si, i — 0,1,2 is a speed depending only on the 
spin of the mode. This is because the aether terms in 
the Lagrangian all contain two derivatives. The scalar 
field action contains terms with fewer than two deriva- 
tives, so in general the spin-0 modes are dispersive. The 
dispersion relation is obtained by taking the determinant 
of Eqs. (24)-(27l considered as a matrix equation. This 



can be simplified to 



3k c 



I 1 lZ Z\( Z Z}Z\ 1 — »L/ -z Z "U Zi Z n 

(uj — k — m )(uj — s Q k ) H —fj, uj /1 k = 0, 

2 C123 

(28) 

where kq is given by (pi) and so is the wave speed of the 
spin-0 mode when /i = [T7] 



(2 - Ci 4 )ci23 



(1 - Ci 3 )(2 + C13 + 3c 2 )ci4 



(29) 



The fact that there are just two, rather than four, prop- 
agating spin-0 modes is a result of diffeomorphism in- 
variance, and is reflected by the fact that Eq. (28) is 
quadratic in uj 2 and therefore admits two solutions. In 
general these modes are superpositions of the spin-0 parts 
of the metric, aether and scalar fields. Their polarizations 
can be obtained by solving Eqs. (24)-(27) with fixed uj 
and k; they depend on the wavenumbcr and are valid 
only in our chosen gauge. 

In the homogeneous limit k — > the dispersion relation 
has the form 



2/2 2,3 2\ n 

UJ (UJ — TO + §KCM ) = 0, 



(30) 



which shows that one mode remains gapless, while the 
mass of the other mode is modified by the \i term. 



IV. LINEARIZED STABILITY AND 
CHERENKOV CONSTRAINT 

In what follows we will consider stability of the classical 
theory linearized about flat space-time. Although this is 
not sufficient to guarantee stability of the full interacting 
quantum theory, or even of the linearized theory about 
a curved background, it is a natural physical condition 
to impose and will lead to constraints on the parameters 
ci.2,3,4, m , M, (i. In the following three sections we will 
show that the linearized theory is stable provided \i is not 
too large compared to to, and the aether parameters lie 
within the allowed parameter range for Einstein- aether 
theory. 

To say a theory is stable is to say that regular "initial 
data" does not grow exponentially in "time" . This notion 
presumes that the theory admits an initial value formu- 
lation with respect to some foliation of the space-time by 
surfaces, called Cauchy surfaces. If the theory is stable 



with respect to one fiducial such foliation, then that suf- 
fices to establish stability from a physical point of view. 
This is because regular initial data on any other Cauchy 
surface can be evolved back to the fiducial Cauchy sur- 
face, where it will again be regular, hence the stability of 
its complete evolution is assured. 

In the case of linearized Einstein-aether theory, the 
field equations are of hyperbolic type, but not with re- 
spect to the space-time metric, since different modes 
travel at different speeds. Nevertheless there is a max- 
imum speed for the modes, and this determines the 
class of time functions whose constant-time surfaces are 
Cauchy surfaces for the full set of equations of motion. In 
particular, the surfaces orthogonal to the aether vector 
in a Minkowski background can serve as Cauchy surfaces 
for the purpose of establishing stability. 

We should note that Carroll et. al. [IH] have proposed 
a different criterion for stability, namely that linearized 
perturbations that are oscillatory on any Lorentz time 
slice should not grow with respect to the corresponding 
time coordinate, and they showed that this criterion is 
strong enough to rule out Einstein-aether theory for most 
values of the parameters c^ .2,3,4- We believe that this 
criterion is not the physically relevant one, because it 
corresponds to imposing stability with respect to a time 
function whose constant-time surfaces are not Cauchy 
surfaces for the equations. This issue will be discussed 
more fully in a forthcoming publication 19]. 

We therefore impose stability in the aether frame. This 



means that the dispersion relation ( 28 ) must have two 
real solutions uj± for every real wavenumber k. Equation 
(28) is a quadratic equation in uj 2 with real coefficients; 



its roots uj 2 . are real provided the discriminant is posi- 
tive. This is necessary but not sufficient, since uj± must 
be positive in order for uj± to be real. We will use the 
fact that two real numbers uj 2 . are both positive if and 
only if the sum uj 2 _ +uj 2 _ and product uj+uj 2 . are both posi- 
tive. This will be convenient to impose since the sum and 
product of roots can be read directly from the coefficients 
of the dispersion relation. In what follows we will work 
backward, first requiring the sum and product roots to be 
positive and then returning to require the discriminant 
to be positive. 

The sum and product of the roots are given by 



2 2 

UJ + UJ_ 



J 2 _ = (l + sl)k 2 + (to 2 - |k cM 2 ) , 
4k 4 + (m 2 - f i 2 /c 123 )s 2 k 2 . 



(31) 
(32) 



These will be non-negative for all k if and only if 



so>0, 


(33) 


KC^JL 2 < 2to 2 /3, 


(34) 


M 2 /ci23 < rn 2 . 


(35) 



Finally, we impose the condition that the roots uj 2 . of 
the quadratic (28) are both real. Thus we require that 



the discriminant be positive, 

(l-s 2 ) 2 fc 4 + 2[(l + s 2 )(™ 2 



,-KcM 2 ) 



2(m 2 - iJ, 2 /c 123 )s 2 ]k 2 + (m 2 - \k cV ?) 2 > 0. (36) 



The discriminant (361 is of the form ak 4 + bk 2 + c with 
a, c > 0. It will be positive if b > or b 2 — 4ac < 0, in 
other words if b > —2\fac. This condition is 



(1 + s\ + \ 1 - s 2 \)(m 2 - l^cn 2 ) - 2{m 2 - M 2 /c 123 )s 2 > 0, 

(37) 



where we have made use of (34) to simplify the expres- 



sion. This inequality can be simplified into two cases 
depending on the sign of 1 — s 2 : 

1/C123 > 3k c /2 
m 2 {l - s 2 ) > (3k c /2 - s 2 /c 123 )^ 2 



if s > 1, (38) 
if s < 1- (39) 



We now consider a further physical constraint on the 
theory: If the phase velocity of scalar aether waves is less 
than the speed of light, then highly energetic particles can 
Chcrcnkov radiate in vacuum. Observations of cosmic 
rays strongly constrain this behavior 20 . We therefore 
require the phase velocity of aether waves to be > 1 for 
all fc. 

First we show that if a mode has phase velocity > 1 for 
some k then its phase velocity remains > 1 for all k. If 
a mode becomes subluminal at a particular wavenumber 
k it would have to satisfy the dispersion relation with 
uj 2 = k 2 , 



r 2/ 2 
[to (s 



l) + (f«c 



s 2 /c 123 )fi 2 }k 2 = 0. 



(40) 



If this equation holds for any k > then it holds for all 
k. This means that if a given mode has speed > 1 for 
any k, it has speed > 1 for all k. 

It is therefore sufficient to enforce that the modes are 
superluminal in the limit k — > oo. To find the phase 



velocities in this limit, we express (28) as a polynomial 
in uj and keep only the dominant power of k in each 
coefficient. The resulting mode speeds are simply 1 and 
So, so the Cherenkov constraint is simply s > 1, as it is 
in pure Einstein-aether theory. 



Thus when the Cherenkov constraint is imposed, (38) 
is the condition for stability. Equations ( 34 ) , ( 35 1 and 



(38) together with Sq > 1 are therefore necessary and 



sufficient conditions for stability and absence of vacuum 
Cherenkov radiation. 

We now derive the condition for the group velocity of 
the linearized waves to be positive. This will play a role in 
the energy positivity constraint. In order for the waves to 
have positive group velocity, we must have doj 2 /dk 2 > 0. 
In the limit k — > oo, the group velocity is positive. Since 
duj 2 /dk 2 is a continuous function of k, it is sufficient to 
show that there is no solution of the dispersion relation 
for which du 2 /dk 2 = 0. 

Solving duj 2 /dk 2 = yields 



2s 2 k 2 = (M7ci23 - m 2 )s 2 + (1 + s z )oj 2 



(41) 



Solving for k, and substituting back into the dispersion 
relation gives aw 4 + bco 2 + c = where 



a = (s 2 - l) 2 , 

b = 2sq[(1 - Sq)to 2 + (1 + So)M 2 / c i23 - 3k c ^ 

1 2 2 / \1 1 

C= (TO - H /Ci23) S . 



(42) 
(43) 
(44) 



Assuming (35) and s 2 > implies c > 0. Just as in the 



discussion following ( 36 1 , there are no real roots for uj 
(i.e. no positive roots for uj 2 ) if and only if b > —2yfac, 
where 



b + 2^fa~c = 4s 2 ,(l/ci23 - 3n c /2)fi 2 



(45) 



Therefore the group velocity is positive for all k if and 
only if Eqs. (33) and (38) hold, as required by stability. 



V. POSITIVITY OF ENERGY DENSITY 

In addition to real frequencies, we further require that 
the linearized perturbations have positive energy. In do- 
ing so we necessarily run into the issue that there is no 
suitable local covariant expression for the energy density 
of a diffeomorphism-invariant theory that accounts for 
the energy in the gravitational field. However for wave- 
like linear perturbations that are periodic in time, it is 
possible to define the average energy density. The pos- 
itive energy conditions for Einstein-aether theory have 
been found both by pseudotensor methods [21] and us- 
ing the Noether current [TB]. The two methods can be 
shown to produce equivalent results [22] ; we will follow 
the latter approach, which is simpler both conceptually 
and computationally. 

The Noether current is defined as follows [22] • We first 
define the canonical one-form 9 a via the total divergence 
that arises when varying the action 



SS = 



g\d 4 x(E&]-5^ + V a 9 a [5iP}) 



(46) 



where E is the equation of motion, and tp schematically 
denotes collectively all the dynamical fields. The Noether 
current 1-form J a [£] associated to the vector field £ a is 
then given by 



j a [$ = e a [Cti>]-L?, 



(47) 



where C^ is the Lie derivative, and L is the Lagrangian. 
When the field equation E = is satisfied the Noether 
current J a constructed this way is conserved, 



V a J a = 0. 



(48) 



The energy is the Noether charge associated with the 
asymptotic time translation t a . It is obtained by choosing 
£ a to coincide with t a at infinity and integrating the flux 
of the corresponding current J a through a Cauchy surface 
E, 



£ = / J a d£ a , 



(49) 



where dE a is the induced volume form on S. 
The Noether current is a sum of two terms Jg- 



and J,, 



and L v re- 



which correspond to contributions from L a 
spectively. The Noether current J x was found by Foster 
[33] . To evaluate J v it is useful to integrate the linearized 
Lagrangian by parts, giving a new Lagrangian 



L'=AV ^VV 



12 2 



nMu a V a f. 



(50) 



With suitable asymptotic boundary conditions this La- 
grangian leads to the same equations of motion. It also 
gives the same Noether charge, and therefore the same 
time-averaged energy density, even though it does not 
give the same Noether current. The advantage of L' over 
L v is that it contains no derivatives of u or g, so the cor- 
responding canonical one-form contains only terms pro- 
portional to Sp, not to Su or Sg. Varying the action, we 
find the contribution of L' to the canonical one-form is 



w^] = (vv + #/)^. 



(51) 



The corresponding Noether current is then determined 
from pfl, 



JZ[S] = (y , <P + nMu a )<p-L'£ a . 



(52) 



We use the Noether current to find the energy density 
of linearized waves following Foster |16j . Consider a com- 
pact source in an asymptotically flat space-time. Fixing 
R to be a sphere of large coordinate radius r, the rate at 
which energy is radiated from the source is given by the 
flux of the Noether charge 



£ = I J-dA, 

IR 



(53) 



where J is the spatial part of J a . The energy is being 
carried away by waves with average energy density u and 
group velocity v g (which we assume positive) , so that the 
average rate of energy loss is 



-(£) = 



uv„ dA, 



(54) 



where (•) denotes time averaging over one period of the 
wave. Equating the two expressions for the rate of energy 
loss we find 



U = (J r ) /Vg 



(55) 



To find the total energy density we need to evaluate this 
formula with the total Noether current, J v + J x . 

To carry out this calculation, note that for sufficiently 
large r all of the dynamical fields can be approximated 
by spherical waves for which 



diip = -{l/v p )4>h, 



(56) 



where v p is the phase velocity of the mode and r 



is the 
Equation (j56h allows spa- 



outward facing unit normal. 

tial derivatives in the energy density to be exchanged for 



time derivatives. The contribution to the energy density 
coming from J v is 



U V = (IP 2 ) /(VpVg). 



(57) 



To find the total energy density we add this to the con- 
tribution to the energy density from the Lagrangian L m 
[Ttj] . Up to a positive factor we obtain 



u ex M 2 (2 - c M )fc 4 if) + 8ci4 (p 2 ) . (58) 

It follows that < C14 < 2 is a sufficient condition for 
positive energy density. We cannot yet conclude that it 
is a necessary condition, since ip and / are related by the 
linearized equations of motion (24127), which imply for 
the complex amplitudes 



[Cl23^ 



n 2 )+v 2 ]ip=hiM(l-c 13 )iujk 2 f. 



(59) 
Substituting this into ( 58 ) we find that up to a positive 



multiplicative factor the energy density is given by 

(2-Cu) [ci2 3 (^ 2 - k 2 - m 2 ) + / i 2 ] 2 + 2c 14Ai 2 (l-c 13 ) 2 w 2 , 

(60) 
where to and k are related by the dispersion relation ( p8| . 
This expression holds for both of the spin-0 modes, each 
of which corresponds to a different solution of the disper- 
sion relation. 

We can now consider two separate limits of the dis- 
persion relation. In the homogeneous limit uj = k = 0, 
the second term in (60) drops out, and positivity requires 
C14 < 2. In the limit k — > oo with uj = k, the first term is 
sub-leading in k, so that C14 > is also needed. The up- 
shot of this calculation is that the constraint from energy 
positivity of the spin-0 modes is 



< C14 < 2, 



(61) 



which is the same as in Einstein-aether theory without 
the scalar coupling. 



VI. COMBINED CONSTRAINTS 

In this section we first gather together all the con- 
straints derived from stability, Cherenkov radiation, and 
energy positivity. Next we consider the constraints from 
post-Newtonian effects, and finally combine these with 
the other constraints to determine the allowed parame- 
ter region. 



A. Stability, Cherenkov and energy constraints 

The constraints discussed already, together with those 
for the spin-1 and spin-2 modes [12J, which are not af- 



fected by the presence of the scalar field, are as follows: 

,2 



constraint ( 62 ) shows that ip cannot be massless if the 



spin-0 stability 
spin-0 speed 
spin-1 speed 
spin-2 speed 



Bkc 1 nv 

< < — S") 

2 C123 H 

(2 - Ci 4 )Ci23 

(1-Ci 3 )(2 + Ci3 + 3C 2 )C14 

2ci -c 2 + c 2 



2ci 4 (l - C13) 

f 

>1, 



>1, 



(62) 

>1, (63) 

(64) 



(65) 

1 -C13 

spin-0 energy < Cu < 2, (66) 

spin-1 energy (2ci - c 2 + cjj)(l - c 13 ) > 0, (67) 



spin-2 energy M > 0. 



(68) 



cosmological Newton constant is to be positive. In order 
to simplify our analysis, we will further assume that m is 
large enough so that the Compton wavelength of <p parti- 
cles is much smaller than any scale on which gravity has 
been tested in the late universe. Newton's law of gravita- 
tion has been tested to sub-millimeter scales, so we shall 
assume the lower bound m > 10 _3 eV on the mass. If 
the mass is smaller than this, a more complete analysis 
would be required to determine the effects of the scalar 
field coupling to gravity. Most likely the mass of the 
inflaton must be tremendously larger than this anyway. 
Assuming that all fields vary on scales much larger 
than 1/m, we can integrate out the scalar field. We first 
solve the linearized scalar field equation ([9]) by an expan- 
sion in D/m 2 , 



Since the spin-1 and spin-2 modes satisfy a dispersion 
relation of the form uj 2 — s 2 k 2 , the stability constraint ip = — (D + m 2 ) -1 [lMQ = — 

for each of these modes is implied by the Cherenkov con- 
straint s 2 > 1. By contrast, the spin-0 modes have a k 
dependent phase velocity, so the separate stability and 
Cherenkov constraints involve several conditions. 



[iM (jlM 



ue- 



, (71) 



The first inequality in ( 62 ) is actually implied by the 
other inequalities. To see this, note that (63 1, (65), and 
(|66|) imply that C123 and kc (m) must have the same sign. 



where in what follows we will keep only the leading or- 
der contribution, ip = —(^M/m 2 )8. If we substitute this 
value back into the linearized scalar field action and ex- 
pand to 0(cp 2 ) we find 



The first inequality in ( 62 ) therefore reduces to the in- 



2m 4 



M 2 M 2 
2m 2 



(72) 



equality C13 < 1, which is already implied by (65) 



B. Post-Newtonian parameters 

As we mentioned in Sec. |IJ Einstein-aether theory has 
a Newtonian limit with a renormalized gravitational con- 
stant Gn given in ^ . The spin-2 energy constraint says 
that G must be positive. Combined with the spin-0 en- 
ergy constraint, this also implies that kjv and therefore 
Newton's constant Gn (p)} is also positive. 

The post-Newtonian parameters /? and 7 take the same 
values (unity) as in general relativity 24J. The effects of 
Lorentz violation are captured at first post-Newtonian 
order by the dimensionless preferred frame parameters 
u\ and a 2 . These can be expressed in terms of 01,2,3,4 as 



Oil 



a 2 = 



The preferred frame parameters are constrained observa- 
tionally to be small [57], ai < 10~ 4 and a 2 < 4 x 10~ 7 , 
and because Einstein-aether theory has four free pa- 
rameters they can be made to vanish exactly in a two- 
dimensional subspace of parameter space. We will see 
that similar considerations apply also in the presence of 
the scalar field. 

We now consider the effect on the post-Newtonian ex- 
pansion of coupling to the scalar field ip. The stability 



The (V0) 2 term is of higher order in the derivative expan- 
sion, and so will be neglected. The remaining term can 
be absorbed in the c 2 term of the aether action, which 
is given by — \M 2 c 2 B 2 . Therefore, upon integrating out 
the scalar field, the new Lagrangian is simply L x with a 
new value of the parameter c 2 which we denote c 2 , 



c-i 



c 2 



/' 



(73) 



The higher derivative terms in the expansion ( 71 ) lead 



to higher derivative terms in the aether action. If we ne- 
glect these higher-order corrections, the post-Newtonian 
parameters can be determined from the known results for 



-8(c§ + 

2c x - c 2 


C1C4) 

+ 4' 










(69) 


Oil 


(ci 


+ 2c 3 


-c0(2c 


+ 3c 2 


+ 


C34) 


(70) 


2 






ci 23 (2- 


cu) 







Einstein-aether theory ( 69 1 and ( 70 ) using c 2 in place of 
c 2 . 

Just as in pure Einstein-aether theory, we can set both 
preferred frame parameters ai, a 2 to zero, now by choos- 
ing 



('2 



Ci 



"2c 2 - C1C3 
3f'i 



' 9 : 



('1 



(74) 
(75) 



Once this choice of C4 has been made, the spin-1 and 
spin-2 speed constraints and the spin-0 and spin-1 en- 
ergy constraints (64)-(67) can be conveniently expressed 



in terms of c± = Ci ± C3, and they hold if and only if 
< c + < 1 and c_ > [55]. Moreover, once this choice 
of c 2 has been made, we have 



C123 



ill 
3ci 



/' 



(76) 



0.6 - 



0.4 



0.2 




As explained above, when the PPN parameters ai % 



c+ 



FIG. 1: The allowed ranges for the parameters 

c± = c\ ± C3 , assuming c 2 and C4 are chosen so that the 

preferred frame PPN parameters vanish. The spin-0 

Cherenkov constraint is satisfied below the solid lines, 

which correspond (from lowest to highest) to /i = 0, 

/1 = |m, \l = \m : fi = m, and /i/m — > 00. The dashed 

line corresponds to the parameter values for which the 

damping rate from weak-field binary systems calculated 

in [16] agrees with general relativity. 



which is positive when the other constraints hold, since 
the positivity of c\ is implied by c± > 0. The second 
inequality of ( 62 1 then becomes 



(i 2 < ci 23 m 2 = (i 2 



^ 13 -m 2 , 



3ci 



(77) 



which is also automatically satisfied when the other con- 
straints are imposed. 

The condition for absence of vacuum Cherenkov ra- 



diation into the spin-0 mode ( 63 1 depends on c 2 , and 
therefore on fj, via ( 74 ) . In pure Einstein-aether theory 
(i.e. when [i = 0), the value of c 2 in Eq. ( 74 ) goes to zero 
as c± —¥ 0. In this limit, sq (29) depends on the ratio 
between c+ and c_. However, if [i 7^ 0, then c 2 does not 
go to zero, and sq — > 00 as c± — > 0. This means that the 
Cherenkov constraint can always be satisfied in an open 
neighborhood of c± — 0. As shown in Fig. Ill the range 
of allowed c± is significantly enlarged as fi increases. 



Other constraints 



have been set to zero by the choices (74) and (75), it 



follows from the other constraints that C123 > 0. Then 
(63) also implies that nc (pi) is positive. This would be 
the value of the cosmological Newton constant relevant 
during inflation. For later times, after the inflaton is in- 
tegrated out 1 , the relevant value is G' c = k' c G, where 
k'q is given by (I6J) but with c 2 in p lace of c 2 . When c 2 
and C4 are set according to ( 74|75 1 it follows [3S] that 
n'c = k n ■ Predictions of primordial element abundances 
from big bang nucleosynthesis are sensitive to G' Cl and 
agreement with observations requires Q3] that G' c not 
differ too much from earthbound measurements of New- 
ton's constant, 



\G' C /G 



N 



II < 1/ 



(78) 



This constraint is thus trivially satisfied when the pre- 
ferred frame PPN parameters vanish. 

The rate of change in orbital period of binary pulsars 
due to gravitational and aether radiation, ignoring ef- 
fects from strong self-gravity of the bodies, was computed 
by Foster [16]. The observed orbital decay agrees with 
that found in general relativity to within better than one 
percent [57]. The constraint that the rate in weak field 
Einstein-aether theory agree with that found in general 
relativity is indicated in Fig. [I] The strong-field effects 
that contribute to the damping rate and orbital dynamics 
of binary systems was computed in |28) in terms of as yet 
unknown parameters that characterize the velocity de- 
pendence of the energy (or action) of the bodies. Strong- 
field effects arc likely to lead to further constraints on the 
parameters, restricting c± to be less than ~ 0.01 — 0.1. 



VII. HOMOGENEOUS COSMOLOGY 

We consider now the cosmology of Einstein-aether the- 
ory coupled to a single scalar inflaton field, with a po- 
tential that depends on the expansion of the aether. The 
metric is assumed to be homogeneous, isotropic, and spa- 
tially flat, hence it can be written in the form 



ds 2 = dt 2 - a(t) 2 dx 2 



(79) 



In a homogeneous and isotropic solution the aether must 
be aligned with the cosmological rest frame. In this 
case, the spatial components of the aether equation are 
automatically satisfied, and the time component of the 
aether equation of motion (jTOj) just determines the La- 
grange multiplier A. That the aether relaxes to such an 
isotropic configuration during cosmological expansion is 



Further constraints arise from the effective cosmolog- 
ical value of Newton's constant in the Friedmann equa- 
tion, and radiation and orbital dynamics of compact bi- 
naries. 



1 The derivative terms in the Klein-Gordon equation are presum- 
ably negligible when H <^C m, which occurs at a temperature 
T <C vrnM. This is satisfied at the nucleosynthesis temperature 
~ 0.1 McV as long asm» 10~ 18 eV. 
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shown perturbatively in [31 [7] . A nonlinear extension of 
this analysis is currently underway j29j . 

We note that our results for homogeneous cosmology 
apply as well to the extended version of Hofava gravity 
[3D] proposed in Ref . [3T] . As shown in Ref. [32] , this fol- 
lows because the aether in these solutions is hypersurfacc- 
orthogonal. 

With the metric ( 79 ) , the expansion is 8 = 3H and the 
homogeneous scalar (Tnflaton) field equation Q is 



tp + 9tp + V v = Q, 



(80) 



where the driving force — V v is now dependent on the 
expansion 9. The terms in V v containing positive powers 
of 9 lead to driving forces that are most relevant during 
inflation when the expansion is large. As the expansion 
slows, these forces subside, allowing for a graceful end to 
inflation. 

The Friedmann equation can be derived by consider- 
ing a metric of the form ds 2 = N(t) 2 dt 2 — a(t) 2 dx 2 and 
varying the action with respect to N. The normalization 
condition for u together with the symmetry completely 
fixes u to u a — (1/N, 0, 0, 0), and the symmetry reduced 
Lagrangian density is (up to a total derivative) 



-gL = Na J 



Ml, , 
— ^9 2 



1 



2N 



tf 



V(9,<p)), (81) 



where = 3H/N and we have introduced a "reduced cos- 
mological Planck mass" Mc = M/^/kc with kc defined 
in (pi). For terms in the action that are linear in 9, the 
factor of 1/N in 9 cancels the factor N in the determi- 
nant of the metric, so that such terms do not contribute 
to the Friedmann equation. This is in agreement with 
the general form of the stress tensor (14). 



Varying (81) with respect to N gives the Friedmann 
equation 



9 2 



M 2 



¥ 2 



V -OVa 



(82) 



where the gauge condition N = 1 has been adopted. This 
equation determines 9 only implicitly as a function of tp 
and ip, since the potential V also depends on 9. 

We now focus on the case when the potential has only 
quadratic terms, 



V(9,<p) 



12 2 



fiM6tp, 



(83) 

) then has the 
e the 9tp cou- 



with pL > 0. The Friedmann equation (82 

standard form for a massive scalar. Whi 

pling does not show up in the Friedmann equation, it 

does affect the Einstein equation, via the pressure term 



in the scalar stress tensor ( 14 ) 



p = V e = fiMtp. 



The infiaton field equation ( 80 1 becomes 
tp + dtp + m 2 ip + jjlMO = 0. 



(84) 
(85) 



The last term acts as an external force that pushes <p in 
the negative direction. 

Slow roll solutions can be obtained by neglecting the 
kinetic term tp 2 in ( 82 ) , and the term ip in the Klein- 
Gordon equation ( 80 1 . The slow roll equations in this 
case become 



3 m 

2Mc 

2 



m 



p 



mMc sgn(p) 



f'c 



(86) 

(87) 



where we have defined 



He 



3kc 



2 + C13 + 3c 2 



This slow roll solution need not be stable: If the Hub- 
ble force —\xM8 is dominant then there is a feedback 
effect, where this force leads to larger field values, which 
leads to a more rapid expansion and therefore a stronger 
Hubble force. The slow roll solution will be stable pro- 
vided that <p and <p are of the opposite sign for all values 
of ip, i.e. if \fj,\ < /i c . This condition is implied by the 
condition for stability of the linearized modes ( 62 ) . This 



is not surprising, since that condition was originally in- 
ferred from the zero wavevector limit of the dispersion 
relation (34), and we have restricted even in the nonlin- 



ear analysis to a quadratic potential. 

The general behavior of this dynamical system for dif- 
ferent values of \i/ /i c is illustrated with phase portraits 
in Fig. [2] Plotted there is the flow on the (tp, <p) plane 
that is obtained when 9 is eliminated using the Fried- 
mann equation. In terms of the dimensionless variables 
tp = tp/Mc and t — mt, we can write the system in the 
form 



where the prime denotes derivative with respect to t. The 
right hand side is the vector field plotted in the figure. 

The number of e-folds of expansion during inflation is 
N = m(a//a,) = J da/a, where Oj and af are the scale 
factors at the beginning and end of inflation, respectively. 
N can be computed in terms of the field values at the 
start and end of inflation as 




N 



1 -[9dt= 1 - 
3./ 3 



9 



dtp. 



(90) 



During slow-roll 9 and tp are given by ( 86 ) and ( 87 1 , and 
tp has a fixed sign, so 



N 



9 2 

pi - pj 



AM 2 ; (1 + Bga(<pi)n/nc) ' 



(91) 



where tpi and tpj are the values of the infiaton field at 
the start and end of inflation, respectively. As [i — > fj, c 
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FIG. 2: Phase portraits of inflaton dynamics with the potential V = ^m 2 ^ 2 + [iMOcp, given by the dynamical 
system (89) with fi > 0. (2a) The result for standard slow-roll inflation. Solutions are attracted toward the slow roll 

solution, which appears as two almost-horizontal lines on the plot. At the end of the slow roll period, a reheating 

phase begins where the field undergoes damped oscillation about the minimum of the potential, which appears as an 

inward spiral in the phase portrait. In (2b) the basic picture remains the same, but the value of (p for the slow roll 

solutions is changed. The duration of the slow roll period is increased if ifii < or decreased if ipi > 0. In (2c), /i is 

at exactly /i c , so any configuration with tp < and tp = is a fixed point, corresponding to a de Sitter space-time. 

When fj, is increased beyond the critical value, as in ( 2d ) , the slow roll solution becomes unstable and ip grows 

without bound. Although only portraits with // > have been shown, the system has combined /i — > —jj, and 

ip — > — p symmetry, allowing the phase portraits with /i < to be obtained from the phase portraits with fi > by a 

reflection through the origin. 
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one of the slow roll solutions approaches a de Sitter solu- 
tion with constant negative <p, corresponding to the limit 
N —¥ oo. These solutions are illustrated in Fig.[2c]on the 
negative <f axis. 

The number of e-folds can be made arbitrarily large 
compared to the /i = case, for fixed values of (fij and 
Mc, by taking fj, — ¥ \i c . However, the ratio /x//i c is con- 
strained by the second inequality in ( |62[ ), which is the 
condition that the spin-0 aether mode not grow exponen- 
tially in the long- wavelength limit. Expressed in terms 
of \i c (88), it yields 



V£_ < 3c i3 + 3c 2 
p* ~ 2 + da + 3c 2 ' 



(92) 



It follows that the limit fi — > fi c is achieved only by C13 — > 
1 or by c 2 -> 00. While the possibility that C13 could 
be close to one has not been ruled out, all the mode 



speeds ( 63 )-( 65 ) diverge in this limit, suggesting that it is 
pathological! - Additionally, it is expected that a strong- 
field analysis of binary systems will constrain c± to be 
smaller than order unity. As for the case where c 2 — >• 00, 
as noted in Sec. |VIB| f can be integrated out so that the 



low-energy theory depends only on c 2 , given by (73). It 



is possible to keep c 2 small even while c 2 becomes large 
by taking ji/m large, though this requires a cancellation 
to occur between c 2 and p 2 /m 2 . 

If the preferred- frame parameters a\ and a 2 are con- 



strained to vanish exactly as in Sec. VI B then as ex- 



plained above, the constraint ( 92 ) is implied by the other 
constraints. In this case, c 2 and therefore [i c are deter- 
mined by fi/m and c± . Then /i//x c can be expressed as 



K 



2m 2 c + + c_(l — c+) 
3/i 2 c + + c_ 



(93) 



This approaches 1 only if fj,/m — > 00, or c + — > 1 and 
c_ — > 00. These are the same conditions as inferred 
above without having constrained ai. 2 to vanish. 

A. Cosmological perturbations 

Finally, we briefly discuss the effect of inflaton-aether 
coupling on the spectrum of perturbations generated dur- 
ing inflation. In single-field slow roll inflation the pri- 
mordial perturbations originate from quantum vacuum 
fluctuations of the coupled inflaton-metric mode. The 
introduction of the aether field leads to an additional 
spin-0 as well as a spin-1 mode, each with their own quan- 
tum fluctuations. If these modes are not directly coupled 
to matter and the background is exactly de Sitter, the 
aether modes decay exponentially [3], though they may 
be sourced during reheating by anisotropic stresses 0]. 
These perturbations may grow in power-law space-times 
for certain values of the parameters 01,2,3,4 [!]• 

Introduction of scalar-aether coupling further modifies 
this scenario. There are two spin-0 coupled inflaton- 
aether-metric modes whose dispersion relation is nontriv- 



ial even in flat space. The generalization of this disper- 
sion relation to curved space-time determines both the 
amplitude of vacuum fluctuations and the time at which 
the modes freeze. Moreover, the modified background 
equation of motion determines how the Hubble parame- 
ter changes with time, and therefore the time at which 
modes of different comoving wavelengths reach the Hub- 
ble radius. 

A priori it is not clear whether the spin-0 modes will 
freeze or decay on superhorizon scales. This behavior is 
determined by how the modes couple to the background 
expansion, and is different for a scalar and the spin-0 part 
of a vector. 



VIII. SUMMARY 

We have considered Einstein-aether theory coupled to 
a scalar field via a potential that depends on the local 
rate of expansion of space in the frame of the aether, 
9 = V a M a . This could be an effective field theory de- 
scription of Lorentz violating UV physics in the vacuum, 
for example, at the scale of a fundamental cutoff. We 
have mostly focused on the lowest order term fiM9ip, 
which leads to breaking of time-reversal invariance, non- 
trivial dispersion, and a cosmological dynamics modified 
by a driving force on the scalar proportional to 9. 

Although we are most interested in this model for the 
potentially observable effects of Lorentz violation on cos- 
mology, we have strived to examine all of the theoretical 
and observational constraints that arise when perturbing 
around a locally flat space-time. We find that when the 
post-Newtonian parameters are matched to general rel- 
ativity, the combined constraints on the remaining free 
couplings ci and C3 imposed by stability, positive energy 
density and absence of vacuum Cherenkov radiation, take 
the same form as in Einstein-aether theory, except for the 
spin-0 Cherenkov constraint, which is relaxed (Fig. [I]). 
There is a single constraint on the new parameter p, the 
second inequality of (62), which is automatically satis- 



fied when the PPN parameters match those of general 
relativity (77). 



An aether field uncoupled to matter does not affect 
the dynamics of homogeneous, isotropic, spatially flat 
cosmology except by a renormalization of Newton's con- 
stant G — > Gc $§■ It turns out that the lowest order 
scalar-aether coupling 9ip contributes isotropic pressure 
and no energy density. Therefore it does not affect the 
Friedmann equation, but the driving force it adds to the 
scalar field equation can act either to oppose or accelerate 
the expansion of the universe in a slow-roll scenario, alter- 
ing the number of e-folds of expansion. An unbounded 
increase in the number of e-folds is possible, although 
this requires carefully chosen values of the couplings in 
order to evade current constraints, and may be ruled out 
by future constraints. This illustrates the importance 
of considering all constraints when building cosmological 
models. 
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A next step would be to derive the spectrum of primor- 
dial perturbations in the model. This would involve gen- 
eralizing the treatment of linearized perturbations from 
flat space to a general inflating background. This differs 
from previous work since the two scalar modes — that of 
the inflaton and that of the aether — are mixed, and be- 
cause the slow-roll dynamics is modified. In particular, 
it would be interesting to determine how the spectral in- 
dex and the tensor-to-scalar ratio depend on the coupling 
constants of the theory. 

While we have focused on the coupling to a scalar infla- 
ton, other types of coupling to the expansion are possible, 
allowing for a wide range of Lorentz-violating phenom- 
ena. For example, any scalar operator could appear in 
the Lagrangian multiplied by 9, allowing effects to be 
switched off (or on) as the expansion of the universe 



slows. In this context, it should be kept in mind that 
in the present universe the largest potentially observable 
values of 9 are not due to the Hubble expansion, but 
rather probably occur outside stellar mass black holes, 
where 9 ~ km~ . 
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